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In a letter [1] Peyrard and Bishop investigated the statistical mechanics of a simple lattice model for 
^he denaturation of the DNA double helix. The model consists of two degrees of freedom u n and v n 
(X£o rr esponding to the transverse displacements of the bases, i.e., displacements along the direction of the 
Qiydrogen bonds, which connect the two bases in a pair. The Hamiltonian for the model is given by 
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H = o m (^n + Vn) + n k { U n ~ Un-lf + {v n ~ ^n-l) 2 + V(u n - V n ) , (1) 
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"^vith the potential between base pairs given by the Morse potential 

O : 2 

V(M n -v n ) = D {exp[-a(u„ - v n )] - l} 2 , 

>-^3vhere D = 0.33eV,a = 1.8A -1 , and the following values for the coupling constant k were used: (a) 
ftx 10- 3 eV/A 2 (b) 3 x lO~ 3 eV/A 2 , (c) 4 x lO^eVyA 2 . In terms of the variables x n = (u n + v n )/y/2, y n = 
.•r4u n — v n )/\/2, the evaluation of the configurational partition function using the transfer integral operator 
ruriethod in the thermodynamic limit reduces to solving a Schrodinger-type equation: 
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1 Ues(y,/3)<Pi(y) = ei<Pi(y) ,P 
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where the temperature-dependent effective potential is given by 
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U eS (y,P) = D{eM-aV2y]-l} + — In . (3) 
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In the thermodynamic limit the ground state dominates with eigen-energy 
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and eigenfunction 

MV) = (^ Q ) 1/2 [f|^p «p(-de-^*) exp[-(d - ±) V2ay] . (5) 

We note that the existence of discrete spectrum is only for values of d = (/3/a)(/cD) 1//2 such that d > |. For 
c? < |, equation (0) has only delocalized states. For the given fc's, we have the corresponding delocalizing 
or melting temperatures: (a)k = 2x lQ- 3 eV/A 2 , T d = 349.868K, (b) k = 3 x 10 -3 eV/A 2 , T d = 428.499if, 
(c) k = 4 x 10~ 3 eV/A 2 ,Td = 494.788-fT. Using the eigenfunction, the authors of [1] evaluated the mean 
stretching (y) = (ipo\y\ipo) an d (y 2 ) = (ipo^lfo) beyond the delocalizing temperatures, where the (y) and 
(y 2 ) in this model are undefined. The figures given by the authors show continuous variations. In addition, 
since ((y — (y)) 2 ) > 0, we have (y 2 ) > (y) 2 . The figures in [1] do not reflect this fact. We have evaluated 
(y) and (y 2 ) and show it in FIG. 1. We note from the eigenfunction (|5[) or from FIG. 1 that the hydrogen 
bonds of between the strands become more stretched for larger temperature or smaller k. 
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